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Abstract 

In the energy range of RHIC and LHC the mechanisms of nuclear suppression of char- 
monia are expected to be strikingly different from what is known for the energy of 
the SPS. One cannot think any more of charmonium produced on a bound nucleon 
which then attenuates as it passes through the rest of the nucleus. The coherence 
length of charmonium production substantially exceeds the nuclear radius in the new 
energy range. Therefore the production amplitudes on different nucleons, rather than 
the cross sections, add up and interfere, i.e. shadowing is at work. So far no theoreti- 
cal tool has been available to calculate nuclear effects for charmonium production in 
this energy regime. We develop a light-cone Green function formalism which incor- 
porates the effects of the coherence of the production amplitudes and of charmonium 
wave function formation, and is the central result of this paper. We found a substantial 
deviation from QCD factorization, namely, gluon shadowing is much stronger for char- 
monium production than it is in DIS. We predict for nuclear effects X2 scaling which 
is violated at lower energies by initial state energy loss which must be also included in 
order to compare with available data. In this paper only the indirect J/^fs originat- 
ing from decay of P-wave charmonia are considered. The calculated Xir-dependence 
of J/ty nuclear suppression is in a good accord with data. We predict a dramatic 
variation of nuclear suppression with xp in pA and a peculiar peak at xf = in AA 
collisions at RHIC. 
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1 Introduction 



Charmonium production off nuclei has drawn much attention during the last two decades 
since the NA3 experiment at CERN has found a steep increase of nuclear suppression 
with rising Feynman xp. This effect has been confirmed later in the same energy range 
@, and at higher energy recently by the most precise experiment E866 at Fermilab ||. No 
unambiguous explanation for these observations has been provided yet. With the advent 
of RHIC new data are expected soon in the unexplored energy range. Lacking a satisfac- 
tory understanding of nuclear effects for charmonium production in proton-nucleus collisions 
it is very difficult to provide a convincing interpretation of data from heavy ion collisions 
experiments @, [f| which are aimed to detect the creation of a quark-gluon plasma using char- 
monium as a sensitive probe. Many of existing analyses rely on an oversimplified dynamics 
of charmonium production which fails to explain even data for pA collisions, in particular 
the observed xp dependence of J/^f suppression. Moreover, sometimes even predictions for 
RHIC employ those simple models. It is the purpose of present paper to demonstrate that 
the dynamics of charmonium suppression strikingly changes between the SPS and RHIC 
energies. We perform full QCD calculations of nuclear effects within framework of the light- 
cone Green function approach aiming to explain observed nuclear effects without adjusting 
any parameters, and to provide realistic predictions for RHIC. 

To avoid a confusion, we should make it clear that we will skip discussion of any mech- 
anisms of charmonium suppression caused by the interaction with the produced comoving 
matter, although it should be an important effect in central heavy ion collisions. Instead, 
we consider suppression which originates from the production process and propagation of 
the cc pair through the nucleus. It serves as a baseline for search for new physics in heavy 
ion collisions. 

The present paper is focused on coherence phenomena which are still a rather small 
correction for charmonium production at the SPS, but whose onset is already observed at 
Fermilab and which are expected to become a dominant effect at the energies of the RHIC 
and LHC. One realizes the importance of the coherence effects treating charmonium pro- 
duction in an intuitive way as a hard cc fluctuation which lose coherence with the projectile 
ensemble of partons via interaction with the target, and is thus liberated. In spite of the 
hardness of the fluctuation, its lifetime in the target rest frame increases with energy and 
eventually exceeds the nucleus size. Apparently, in this case the cc pair is freed via inter- 
action with the whole nucleus, rather than with an individual bound nucleon as it happens 
at low energies. Correspondingly, nuclear effects become stronger at high energies since the 
fluctuation propagates through the whole nucleus, and different nucleons compete with each 
other in freeing the cc. In terms of the conventional Glauber approach it leads to shadowing. 
In terms of the parton model it is analogous to shadowing of c-quarks in the nuclear struc- 
ture function. It turns out (see Sect. |j) that the fluctuations containing gluons in addition 
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to the cc pair are subject to especially strong shadowing. Since at high energies the weight 
of such fluctuations rises, as well the the fluctuation lifetime, it becomes the main source of 
nuclear suppression of open and hidden charm at high energies, in particular at RHIC. In 
terms of the parton model shadowing for such fluctuations containing gluons correspond to 
gluon shadowing. 

The parton model interpretation of charmonium production contains no explicit coher- 
ence effects, but they are hidden in the gluon distribution function of the nucleus which 
is supposed to be subject to QCD factorization. There are, however, a few pitfalls on this 
way. First of all, factorization is exact only in the limit of a very hard scale. That means 
that one should neglect the effects of the order of the inverse c-quark mass, in particular 
the transverse cc separation (r^) ~ f/ml- However, shadowing and absorption of cc fluc- 
tuations is a source of a strong suppression which is nearly factor of 0.5 for heavy nuclei 
(see Fig. |). QCD factorization misses this effect. Second of all, according to factorization 
gluon shadowing is supposed to be universal, i.e. one can borrow it from another process 
(although we still have no experimental information about gluons shadowing, it only can 
be calculated) and use to predict nuclear suppression of open or hidden charm. Again, fac- 
torization turns out to be dramatically violated at the scale of charm and gluon shadowing 
for charmonium production is much stronger than it is for open charm or deep-inelastic 
scattering (DIS) (compare gluon shadowing exposed in Fig. |7| with one calculated in || for 
DIS). All these important, sometimes dominant effects are missed by QCD factorization. 
This fact once again emphasizes the advantage of the light-cone dipole approach which does 
reproduce QCD factorization in situations where that is expected to be at work, and it is 
also able to calculate the deviations from factorization in a parameter free way. 

Unfortunately, none of the existing models for J /iff or iff' production in NN collisions 
is fully successful in describing all the features observed experimentally. In particular, the 
J/\&, iff' and Xi production cross sections in NN collisions come out too small by at least an 
order of magnitude 0. Only data for production of Xi whose mechanism is rather simple 
seems to be in good accord with the theoretical expectation based on the color singlet 
mechanism (CSM) [8|, |9| treating xi production via glue-glue fusion. The contribution of 
the color-octet mechanism is an order of magnitude less that of CSM |§, and is even more 



suppressed according to |]ID |. The simplicity of the production mechanism of \2 suggests 



to use this process as a basis for the study of nuclear effects. Besides, about 40% of the 
J/iffs have their origin in x decays. We drop the subscript of X2 in what follows unless it is 
important. 



1.1 What has been understood already 

A lot of work has been done already and considerable progress has been achieved in the 
understanding of many phenomena related to the dynamics of the charmonium production 
and nuclear suppression. We would like to start with reviewing some of these phenomena 
which are employed in present paper. 

• Relative nuclear suppression of J /iff and iff' has attracted much attention. The iff' 
has twice as large radius as the J /iff, therefore should attenuate in nuclear matter much 
stronger. However, formation of the wave function of the charmonia takes time, one can- 
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not instantaneously distinguish between these two levels. This time interval or so called 
formation time (length) is enlarged at high energy by Lorentz time dilation, 



2 Em 

and may become comparable to or even longer than the nuclear radius. In this case neither 



J/ty, nor propagates through the nuclear medium, but a pre-formed cc wave packet [I I 
Intuitively, one might even expect an universal nuclear suppression, indeed supported by 
data 0, |, H]. However, a deeper insight shows that such a point of view is oversimplified, 
namely, the mean transverse size of the cc wave packet propagating through the nucleus 
varies depending on the wave function of the final meson which the cc is projected to. In 



particular, the nodal structure of the 25 state substantially enhances the yield of \J/' |13|, [14 
(see in |15| , |i~6|| a complementary interpretation in the hadronic basis). 

• Next phenomenon is related to the so called coherence time. Production of a heavy 
cc is associated with a longitudinal momentum transfer q c which decreases with energy. 
Therefore the production amplitudes on different nucleons add up coherently and interfere 
if the production points are within the interval l c = l/q c called coherence length or time, 

t -l^L ( 2 ) 

This time interval is much shorter than the formation time Eq, ([!]). One can also interpret it 
in terms of the uncertainty principle as the mean lifetime of a cc fluctuation. If the coherence 
time is long compared to the nuclear radius, t c > Ra, different nucleons compete with 
each other in producing the charmonium. Therefore, the amplitudes interfere destructively 



leading to an additional suppression called shadowing. Predicted in [13] this effect was 



confirmed by the NMC measurements of exclusive J/^f photoproduction off nuclei [O (see 



also |14|])- The recent precise data from the HERMES experiment [|T^] for electroproduction 



of p mesons also confirms the strong effect of coherence time [19| 



Note that the coherence time Eq. (fj) is relevant only for the lightest fluctuations |cc). 
Heavier ones which contain additional gluons have shorter lifetime. However, at high ener- 
gies they are also at work and become an important source of an extra suppression (see in 
and Sect. §). They correspond to shadowing of gluons in terms of parton model. In terms of 
the dual parton model the higher Fock states contain additional qq pairs instead of gluons. 
Their contribution enhances on a nuclear target lead to softening of the xf distribution of 
the produced charmonium. This mechanism has been used in [^U] to explain xf dependence 
of charmonium suppression. However the approach was phenomenological and data were 
fitted. 

The first attempt to implement the coherence time effects into the dynamics of char- 



monium production off nuclei has been done in [21]| . However, the approach still was phe- 
nomenological and data also were fitted. Besides, gluon shadowing (see Sect. |]) has been 
missed. 

• The total J/\l/-nucleon cross section steeply rises with energy, approximately as s a2 . 
This behavior is suggested by the observation of a steep energy dependence of the cross 
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section of J/ty photoproduction at HERA. This fact goes well along with observation of 
the strong correlation between XBj dependence of the proton structure function at small 
Xsj and the photon virtuality Q 2 : the larger Q 2 is (the smaller is its qq fluctuation), the 
steeper the F 2 (xBj,Q 2 ) rises with 1/xsj- Apparently, the cross section of a small size 
charmonium must rise with energy faster than what is known for light hadrons. The J/ty- 



nucleon cross section has been calculated recently in |22j employing the light-cone dipole 
phenomenology, realistic charmonium wave functions and phenomenological dipole cross 
section fitted to data for F 2 (x, Q 2 ) from HERA. The results are in a good accord with data 
for the electroproduction cross sections of J/ty and ty' and also confirm the steep energy 
dependence of the charmonium-nucleon cross sections. Knowledge of these cross sections 
is very important for understanding of nuclear effects in the production of charmonia. A 



new important observation made in p2| is a strong effect of spin rotation associated with 
boosting the cc system from its rest frame to the light cone. It substantially increases the 
J /\& and especially photoproduction cross sections. The effect of spin rotation is also 
implemented in our calculations below and it is crucial for restoration of the Landau- Yang 
theorem (see Appendix C). 

• Initial state energy loss by partons traveling through the nucleus affects the xp distribu- 
tion of produced charmonia [23] especially at medium high energies. A shift in the effective 



value of Xi, which is the fraction of the incident momentum carried by the produced char- 
monium, and the steep Xi-dependence of the cross section of charmonium production off a 
nucleon lead to a dramatic nuclear suppression at large x\ (or xp) in a good agreement with 
data [[I], 0]. The recent analyses J24], [25| of data from the E772 experiment for Drell-Yan 



process on nuclei reveals for the first time a nonzero and rather large energy loss. 
1.2 Outline of the paper 

The paper is organized as follows. The light-cone (LC) formalism for gluon-nucleon collision 
formulated in the rest frame of the target nucleon is introduced in Sect. |[ As usual, the 
amplitude is represented as a convolution in the impact parameter space of the LC wave 
functions of the incident gluon and the final charmonium, and the dipole transition cross 
section. The latter corresponds to the interaction with the target which causes a transition 
between \cc)g and \cc) x , which are the charm-anticharm pairs with quantum numbers of a 
gluon and a x respectively. The transition amplitude S tr can be expressed in terms of the 
ordinary flavor independent dipole cross section aq q of interaction of a colorless qq dipole 
with a nucleon which is rather well known from phenomenology. 

While the perturbative gluon wave function, an analog to the photon one, is known, the 
LC wave function of a charmonium needs to be constructed. Even if the wave function in 
the rest frame of the charmonium is known, it is not a straightforward procedure to boost it 
to the infinite momentum frame. We apply the widely accepted prescription for the Lorentz 
boost, which is rather successful in describing data for exclusive photoproduction of char- 
monia [ 22 1 . It is presented in Appendix A for the case of the \ wave function. One may 
not expect important relativistic effects for such a heavy system as a charmonium. Indeed, 
it is demonstrated in Appendix B that the distribution over fraction of the x momentum 
carried by the quarks peaks at 1/2 with a tiny width. However, the Melosh spin rotation 
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effect is still very important. It is demonstrated in \22\ how much it affects the photopro- 
duction cross section of J/^> and especially \E f/ . In the case of \ gluoproduction inclusion 
of the spin rotation restores the Landau- Yang theorem which forbids production of by 
an on-mass-shell gluon. Otherwise, it would be badly violated even in the nonrelativistic 
limit where one might expect the spin rotation effects to vanish. What is rather obvious in 
the parton model, still needs special efforts to be proven within the LC approach. The fine 
tuning leading to cancelation of different parts of the \i production amplitude is demon- 
strated in Appendix C. This can be considered as another strong support for the Lorentz 
boost procedure we are using. 

The production of charmonia off nuclei is controlled by few length scales as is discussed 
above and in Sect. Once a colorless cc wave packet is produced, it is evolving during 
propagation through the nucleus due to absorption and transverse motion of the quarks. 
The result should be projected to the charmonium wave function. The length scale of the 
evolution is controlled by the formation time Eq. ([!]) . 

Another much shorter length scale is controlled by the coherence time Eq. (0). At high 
energy it becomes long leading to shadowing which enhances the suppression of charmonia. 
This is the new phenomenon which is not yet at work at the SPS, but is expected to nearly 
saturate at RHIC. 



We start with the effects of quantum coherence for simple limiting situations. In Sect. qTL 



we study the case of a very long If ^> l c ^> Ra when the transverse quark separation is 
frozen by the Lorentz time dilation for the time of propagation through the nucleus. While 
the final state absorption occurs with the conventional dipole cross section cx^ty), it is 



not obvious which cross section controls shadowing. It is the key observation of Sect. 3.1 



that this is the dipole cross section of interaction of a colorless system consisting of three 
partons, \qqG). The combined effect of shadowing and absorption is given by Eq. (^) - 



The case of a not too long coherence time l c ~ Ra is considered in Sect. |3T2| . Since 
If l c , the time taken by the produced cc for its further development can be rather long, 
// ^> Ra and the transverse size of the cc pair can be treated as frozen. Eq. ( f4"5|) interpolating 
between the regimes of very short and long coherence lengths is governed by the longitudinal 
momentum transfer q c = 1/ la- 
in the general case considered in Sect. [O] any of l c and If can be either short or long 
and one must take care of the transverse size fluctuations of a color dipole propagating 
through the nucleus. The appropriate approach is the path integral technique summing up 
all possible paths of the partons Evolution of a cc wave packet in nuclear medium is 
described by the LC Green functions satisfying the two dimensional Schrodinger equations 
(|52|), which are different for colorless and colored dipoles. The central result of this paper, 
the amplitude of the process GA — > \X is described by Eq. ([49]) which is illustrated 
pictorially in Fig. |3|. 

Evaluation of the effects of coherence and formation is performed in Sect. |3.4| and the 
results are demonstrated in Fig. |j. The nuclear transparency first rises with energy due to 
the formation length effects (color transparency), but then, falls down at higher energies 
due to growing coherence length. 

In Sect. H a third scale governing the nuclear effects is introduced. It is related to the 
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higher Fock components in the projectile gluon, \ccG), \cc2G) etc. Since the gluons usually 
carry a small fraction of the total momentum, these fluctuations are rather heavy compared 
to \cc), hence they have a coherence time shorter than in Eq. (|2]). Therefore, the contribution 
of these fluctuation to nuclear shadowing which must be associated with gluon shadowing, 
is delayed down to smaller values of x-i- 

First of all one must develop an impact parameter approach for fluctuations containing 
gluons. It is quite a difficult task, but it pays off when one needs to calculate shadowing 
which is given by a simple eikonalization, or by employing the Green function technique 
if 22 is not small enough. This is done in Appendix D and the results are summarized in 
Sect. |4j]. It turns out that the twelve Feynman graphs for the LC wave function of the 
\ccG) depicted in Fig. ^ and calculated in Appendix D are reduced to the main contribution 
which corresponds to the production of x directly from the projectile gluon prior or after 
the interaction with the target, as is illustrated in Fig. || Therefore the dipole formalism 
appropriate for gluon shadowing is pretty similar (up to color factors) to that for Drell- 
Yan reaction |26|, |2?|, |2£| 59|. The LC wave function of the \ccG) state where the cc pair 



is colorless, turns out to be very different from that for the \qqG) fluctuations in DIS 
where the strong nonperturbative interaction between the color-octet qq and the gluon 
substantially diminishes the qq - glue separation and reduces shadowing of gluons. Indeed, 
gluon shadowing for x production depicted in Fig. |4.2| is much stronger than calculated for 
DIS in 0. 

There are some indications that gluons in nuclei may be enhanced at at X2 ~ 0.1. This 
effect increasing the production rate of charmonia at smallx^ is discussed and estimated in 
Sect. [Dj 

More effects are to be included in order to compare with data. In Sect. |5| the corrections 
for energy loss by the projectile partons which are due to initial state interactions are 
evaluated according to the prescription of [|3|, |4], p5| . This effect dramatically violates 
X2 scaling. While it is the main mechanism shaping the xp dependence of charmonium 
suppression at medium high energies (the SPS energy range), its influence is essentially 
reduced at Fermilab, and no visible energy loss corrections are expected at RHIC. 

Eventually, in order to compare with available data for J/\l/ production off nuclei one 
has to make corrections for decay x J/^ 1- This is done in Sect. ^ and the corrections 
are found to be rather small. Indeed, they are basically the same in pN and pA collisions 
and essentially cancel out in the ratio. 

Incorporating all these effects in our calculations in Sect. |7| we arrive at quite a good 
agreement with available data. While energy loss is the dominant effect at 200 GeV (Fig. |T0|) , 
the steep xp dependence of nuclear suppression at 800 GeV (Fig. |) is a combined effect of 
quark and gluon shadowing and energy loss. This is the first manifestation of the onset of 
gluon shadowing, but it is expected to to become the main source of charmonium suppression 
at RHIC. Our predictions depicted in Fig. |ll] demonstrate a dramatic variation of nuclear 
suppression versus small xp for pA collisions. Translated into suppression in AA collisions 
the same effects form a peculiar narrow peak at xp = 0. 

Our conclusions and outlook are presented in Sect. [| This is the first full QCD calcula- 
tion of coherence effects in charmonium production off nuclei demonstrating that coherence 
is the dominant phenomenon which governs nuclear suppression of charmonia at the energies 
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of RHIC and LHC. 



2 The light-cone dipole formalism for charmonium pro- 
duction off a nucleon 

The important advantage of the light-cone (LC) dipole approach is its simplicity in the 
calculations of nuclear effects. It has been suggested two decades ago PDA that quark con- 



figurations (dipoles) with fixed transverse separations are the eigenstates of interaction in 
QCD. Therefore the amplitude of interaction with a nucleon is subject to eikonalization 
in the case of a nuclear target. In this way one effectively sums the Gribov's inelastic 
corrections in all orders. 

Assuming that the produced cc pair is sufficiently small so that multigluon vertices can 
be neglected we can write the cross section for G N — > %X) as (see Fig. []]), 

a(GN -> X X) = l_ £ j ^a s (k 2 T )Hx,k 2 T ) \M ab (k T )\ 2 , (3) 

\ c / a,b T 



where A;^) = dG(x, k^)/d(ln k^) is the unintegrated gluon density, G(x,kj,) = xg(x,kj> 

'i 



x = M%/s); M a b(kx) is the fusion amplitude GG — > x with a, 6 being the gluonic indexes. 
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Figure 1: Perturbative QCD mechanism of production of the \ states in 
a gluon-nucleon collision. 



In the rest frame of the nucleon the amplitude can be represented in terms of the cc LC 
wave functions of the projectile gluon and ejectile charmonium, 



ab 



M ab (k T ) = -j= 



da / d 2 r T £ (r T , 



a 



$£f(f T ,a) 



(4) 



tie 



where 



n = (1 - a) f T , 



r 2 = —atT ■ 



(5) 
and 

(cc, {1}| from the LC wave function of the gluon and charmonium respectively, and calculate 



For the sake of simplicity we separate the normalization color factors (cc, {8} a 



S 



the matrix element, 

(5 C| {8} |iA 6 |g C ,{l}) = ^ I (6) 

which is shown explicitly in (Q). Thus, the functions ^Q^ x Jfx,a) in (|j) represent only the 
spin- and coordinate dependent parts of the corresponding full wave functions. 
The gluon wave function differs only by a factor from the photon one, 



m?T,a) = ^ &fd^K (er T ) , (7) 
where is the c-quark spinor, and 



f = i* v $ , (8) 
6 = m c a-e + i(l-2a)(a-n)(e-V) + (exn)V , (9) 
e 2 = Q 2 a(l-a)+m 2 c , (10) 



2 -o2 /i \ i 2 



The gluon has virtuality Q 2 and polarization vector e and is moving along the unit vector 
n (in what follows we consider only transversely polarized gluons, e ■ n = 0). 

The expression for the LC wave function of a charmonium related by the Lorentz boost 
from the charmonium rest frame is rather complicated and is moved to Appendix A. This 
complexity is a consequence of the nonlocal relation between the LC variables {fr, ot) and 
the components of the 3-dimensional relative cc radius-vector r in the rest frame of the 
charmonium. Also the Melosh spin rotation leads to a nontrivial relations between the two 
wave functions (see e.g. in f22| ). This is a relativistic effect, it vanishes in the limit of small 
velocity v — *■ of the quarks in the charmonium rest frame. 

A word of caution is in order. In some cases the Melosh spin rotation is important even 
in the limit of vanishing quark velocity v — *■ 0. An example is the Landau- Yang theorem 
||31f which forbids production of the \i state by two massless gluons. However, the LC 
approach leads to creation of the xi even in the limit v — > if the effect of spin rotation 
is neglected. It is demonstrated in Appendix B that the Landau- Yang theorem is restored 
only if the Melosh spin rotation is included. Such a cancelation of large values is a kind of 
fine tuning and is a good support for the procedure of Lorentz boosting which we apply to 
the charmonium wave functions. 

Since the gluon LC wave function smoothly depends on a while the charmonium wave 
function peaks at a = 1/2 with a tiny width estimated in Appendix B, {(a — 1/2) 2 ) = 0.01, 
we can replace the charmonium wave function in the matrix element in (fj) with 

$f(r T ,a)«<5(a-^) J da^(r T ,a) . (11) 

It is convenient to expand the LC charmonium wave function in powers of v. The result 
depends on the total momentum J and its projection J z on the direction n. The charmonium 
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wave function integrated over a has the form 

'da&» (r T ,a) = 
where 



a ■ e± + — (e± x n) ■ V - — ^ (e± ■ V) (a ■ V) 
m 2 



r T 



R(r T ) + -r^R"(r T ) + 0(v 4 ) 
Ami 



ew + o( V 4 ) , 

(12) 
(13) 



and R(r) is the radial part of the P-wave charmonium in its rest frame (see derivation in 
Appendix A). The new notations for the polarization vectors are, 



+ ie y 



(14) 



In what follows we use the LC wave functions of gluons and charmonium in order to 
calculate matrix elements of operators which depend only on the LC variables and a. 
Therefore, for the sake of simplicity we can drop off the indexes fi, ft and summation over 
them, i.e. replace 



£ ($f (f T , a))* ^{r T , a) => $* (f T , a) $ G (f r , a) 



(15) 



With this convention we can rewrite the cross section Eq. (HI) as, 



x 



cr(GiV ^xX) = J da J da' J d 2 r T d 2 r' T 
o o 

{$* (f T , a) $ x (4, a') S fr (f T , 4, a, a') $ G (f T , a) $^(4, a')} , 



where the transition cross section £*' r is a combination of dipole cross sections, 

E' r (f T , 4, a, a') = — [a 5 g(fi - f£) + cr ? ,(f 2 - r[) - cr^n - r/) - o"g g (f 2 - f^) 
lb L 



and fi, and are defined like in Eq. (|5|). The dipole cross section [|3D 

47r /■ d 2 k T 



Vqq{rT, s) 



J^ra s (k 2 T )F(x,k 2 T ) (l- e *"* 



(16) 



(17) 



18) 



corresponds to interaction of a colorless gg pair of transverse separation tt with a nucleon 
at the squared cm. energy s and x = 4k 2 -,/ s. The explicit s-dependence of the cross sections 
is dropped in Eqs. (|i~6|), (|i~7l). 

Since small distances ~ ~ 1/ m c dominate in the integral in Eq. (0) one can make 
use of the approximation 



r T ^0 
10 



C(s) r 2 T , 



(19) 



then S tr reduces to the very simple form, 



E^(f T , 4, a, a') = ^ (f T • 4) = £ (g* • d) (e* ■ d')* , (20) 



A=l,2 



where 



I( m Pt , d ^m^. (21) 



The vectors e t A can be interpreted as a polarization vector of the Weizsacker- Williams gluon 
of the target. 

Within this approximation the cross section ( |16D can be represented in the form, 

a(GN^ X X)=Y / \A W \ 2 , (22) 

A 

where 

i 

A (A) = J da J d 2 r T $* (f T , a) (4 • d) $ G (r T , a) . (23) 



3 Production of charmonia off nuclei: shadowing of 
c- quarks 

Nuclear effects in the production of a \ are controlled by the coherence and formation lengths 
which are defined in ([[]), (|2|). One can identify two limiting cases. The first one corresponds 
to the situation where both l c and If are shorter that the mean spacing between bound 
nucleons. In this case one can treat the process classically, the charmonium is produced on 
one nucleon inside the nucleus and attenuates exponentially with an absorptive cross section 
which is the inelastic x ~ N one. This simplest case is described in [f23], [32) . 

In the limit of a very long coherence length l c 3> Ra one can think about a cc fluctuation 
which emerges inside the incident hadron long before the interaction with the nucleus. 
Different bound nucleons compete and shadow each other in the process of liberation of this 
fluctuation. This causes an additional attenuation in addition to inelastic collisions of the 
produced color-singlet cc pair on its way out of the nucleus. Since l c <C If an intermediate 
case is also possible where l c is shorter than the mean internucleon separation, while If is of 
the order or longer than the nuclear radius. 



3.1 The high energy limit, l c ^> Ra 

We start with the case of very long coherence and formation lengths, l c , If 3> Ra- In this 
case the incident gluon converts into a color-octet cc fluctuation with a lifetime much longer 
than the nuclear radius. Therefore one can consider this fluctuation as produced long in 
advance and propagating through the whole nucleus as is illustrated in Fig. 0. Then we can 
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Figure 2: The incident gluon converts onto a cc pair long in advance 
of the nucleus. The pair propagates and attenuates with the absorption 
cross section 03 (tt) (see the text) up to the point where it converts 
into a colorless cc pair with quantum numbers of X- Then it continues 
propagating through the nucleus, and is attenuated with the cross section 



employ the results of Ref. |H| for the evolution of the density matrix a qq (cc in our case) 
wave packet propagating through the nuclear medium, 



R{f T ,f T ,a,a'\z) = R m (? T ,i%.,a,a/\z)Pi + - R {8) {f T ,f T ,a,a'\z) 



8 , 



(24) 



where R(rr, ?ti a i oc'\z) is the density matrix of the cc propagating through nuclear medium 
which depends on the transverse separation ty, fraction a of the LC momentum, and on the 
longitudinal coordinate z; P\ and P 8 are the projection operators on the singlet and octet 
states of the qq, respectively, and satisfy the relations, 



Pi + Ps = 1 ; 

Tr A = 1 ; Tr P 8 = 8 . 
The matrices and R^ are the solutions of the evolution equations, 



dR^(f T , r^, a, a'\z) 
dz 



Pa{z) 



■-E^R^\f T ,r^,a,a'\z) 



+ R {8 \r T ,r T ,a,a'\z) 



(25) 



(26) 



dR( 8 \r T ,r T ,a,a'\z) 
dz 



E (8) R^{r T ,rlr,a,a'\z) 



a, a \z) 



rp, IX, 



(27) 



where the nuclear density Pa(z) depends on the longitudinal coordinate z and (implicitly) 
the impact parameter b. The transition cross section S ir is defined in (|i~7D . 
The function 

£« = E«(r T ,4,a,c/) = a m {r T ) + a qq (f T ) (28) 
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is the total cross section for the interaction of a 4-quark ensemble, two color-singlet cc pairs, 
with a nucleon. Correspondingly, 



E (8) 



£< 8 >i 



r T ,rfp,a,a') 



1 r 



2 Oqq{n - r£) + 2 cr 5g (f 2 - f-[) 



+ 7 o- ?g (fi - r*i) + 7 a 5g (f 2 - 7*) - <r gg (fi - f 2 ) - a, 



99 V 1 



(29) 



is the total cross section, where a 4-quark system consisting of two color-octet cc pairs whose 
centers of gravity coincide, interact with a nucleon. 

The initial conditions for a cc pair originating from the projectile gluon are, 



r T , r^, a, a'|z) 



rr, r T , a, a \z l 







$ G (f T ,a)$^(4,a') . 



(30) 
(31) 



With the evolution equations (p6j), (|27| ) we are in position to extend Eq. (|T^) to the case 
of a nuclear target, 



1 1 

a(GA -> x-X") = J d 2 b J da J da' J d 2 r T d 2 r' T $*(f T , a) $ x (4, a') 



i?«( 



r T ,r T ,a,a |2; + J , 



(32) 



where 2+ — > 00, and i?^ implicitly depends on the impact parameter b. This expression 
includes all the multiple color exchanges between the projectile cc pair and bound nucleons 
in the target which eventually convert the initial color octet cc into the final colorless state 
we are interested in. 

Color exchanges on different nucleons add up incoherently, and the cross section is rather 
small due to smallness of the mean transverse separation for the heavy cc pair. Therefore, 
we keep only the lowest (first) order in S* r , but all higher orders in S 8 and S 1 . Then, the 
cross section of the process G A ^> \ X reads, 



where 



o(GA -> X X) = j ' d 2b j da J da ' J d 2 r T d 2 r' T 



{%(f T , a) $ x (4, a') S(r r , f* T , a, a') $ G (r T , a) $ G (r^, a')} , 



£(f T , f£, a, a') 



2E^(f T ,f;«,a') 



£« (f T , r^r, a, a') - £( 8 ) (f T , 4, a, a') 



exp 



zW(r T ,r^,a,a')T(b) 



— cxp 



--EW(f T) f>,a')T(6) 



Here 



T(b) 



dzp A (b,z) 



(33) 



(34) 



(35) 
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is the nuclear thickness function. 

Eq. (|34| ) can be modified to a form which makes its physical meaning more transparent 
(for the sake of brevity we drop the variables r T , r' T , a, a' in functions E* r , and E( 8 )), 



E(f T , r' T , a, a') 



d 2 b / dz p A (b, z) exp 



x 



E (tr) 



exp 



(36) 



where 



DC 

T+(b,z) = Jdz'p A (b,z') 

z 
z 

T-(b,z)= J dz' PA (b,z>) 



Furthermore, we introduce the following notations, 



a 3 [r T ,a 



Tt, r T , a, a 
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r' T =0 



- [a qq (ar T ) + a 5(? [(l - a)r T ] 



AE^( 



r T , r^, a, a') 



S (8) ( 



r T , fy, a, a') - <x 3 (f T , a) - a 3 (r^, a') 



(37) 



(38) 
(39) 



Since heavy quarks are expected to be produced with a small separation, we can employ 
the approximation (|K|). 



AE^l 



7Y, r T , a, a 



7 -9(a + a') + 18 a a' 



7Y • r T 



(40) 



This amplitude describes the diffractive transitions of the cc octet accompanied with change 
of the orbital momentum and parity of the cc pair. Since the color-exchange amplitude also 
has this property the product of these two amplitudes corresponds to creation of direct J/ty 
(&'■■■) on a nuclear target without radiation of any extra gluons. In this paper, however, 
we concentrate on the production of x s an d neglect this effect, i.e. replace 



E^i 



tt, r T , a, a 



a 3 (f T ,a) + <7 3 (7Y,c/) . 



(41) 



If one uses the small-r^ approximation (|T^) also for £(* r ), then the cross section Eq. (1.19) 
can be represented in the form, 



a{GA^ X X) = J d 2 b J dzp A {b,z)^2\ AiX) (b, 



(42) 



where 



A™(b,z) = j 'da J d 2 r T <S>* x (r T ,a)A( x \b,z-,r T ,a)<S> G (r T , 



(43) 
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and 



A (x \b,z;f T ,a) 



exp 



d exp 



--«7 3 (r T ,a)T_(6, z) 



(44) 



The vector d is defined in 

In Eq. fl44|) the first exponential factor corresponds to the nuclear attenuation of the x 
produced in the color-exchange rescattering of the projectile cc. The second exponential 
should be interpreted as nuclear shadowing for the diffractive transition G — > cc which 
preserves the initial quantum numbers and color (see in ||34|1 ). The factor • (/corresponds 
to the amplitude of dipole radiation (or absorption) of a gluon by the color-octet cc pair at 
the point z of the color-exchange interaction. 



3.2 Medium high energies, l c ~ Ra, but If ^> Ra 

Such an energy range is possible only for heavy flavor production due to the relation If ^> l c 



| i~3| . In this case we can still treat the final state colorless cc pair propagating through 
the nucleus as having a frozen transverse separation, while corrections for finiteness of l c 
must be done. Correspondingly, the operator A^(b, z; f T , a) has to be modified and can be 
represented as a sum of two terms, 



where 



AM (6, z; Ft, a) = 4 A) (&, *; r r , a) + 4 A) (6, z; f T , a) 



(A)/ 



(fo, z; r T , a) = exp 



iqiz- -aq Q (r T )T + (b, z) 



e< A ) ■ d 



(45) 



(46) 



x 



AY\b, z; f T , a) = -- exp 



iq L z- -ag q (r T )T + (b, z) 



J dz\ exp 



iqizi- -a 3 {r T a)T_(b, z x ] 



ei A) ■ d 



cr 3 (r T ,a) p A (b, z x ) 



(47) 



where qi = 1/Z C is the longitudinal momentum transfer to the nucleus. This expression 
interpolates between the high-energy limit l c 3> Ra {ql — 0) where it acquires the form 
and the limit of l c shorter than the mean nucleon separation when — > 0. Such 



was first derived in p5f for inelastic diffractive 



a form of the amplitude Eqs. (f45|) - 
photoproduction of J/ty off nuclei. 

The term A\ in Eq. d45|) describes the production of the colorless cc pair at (b, z) by 



the incident gluon. This pair then attenuates escaping the nucleus, while the gluon has no 
initial state interaction (leading to production of a cc pair). 

The term A 2 in Eq. fl4"5| ) describes the diffractive production of a cc pair by the gluon 
with the same quantum numbers (except the color dipole moment of the cc) at the point 
Zi < z. This color-octet cc pair propagates and attenuates between the points z\ and z with 
the cross section 03 (ry, a). Then at point z it experiences a color exchange interaction and 
produces a colorless cc pair with the quantum numbers of a %• 
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In fact, Eq. ( fl"5| ) breaks down towards the low-energy limit at x energies of a few tens 
of GeV because If becomes comparable with Ra and one cannot neglect any more the 
fluctuations of the transverse size of the cc pair during its propagation through the nucleus. 
In this case one should apply the technique of the light-cone Green function describing the 
propagation of qq pairs |36|, |37|, |28], ||. 



3.3 General case 

The transition between the limits of very short and very long coherence lengths is performed 
using the prescription suggested in [33| for inelastic photoproduction of J/\l/ off nuclei. The 
cross section of x production off a nucleus can still be represented in the form (f42"f ), but the 
amplitude is modified as, 

1 

A (X \b,z) = J da J d 2 r T J d 2 r' T $* (f T , a) £> {A) (f T , r' Tl a; 6, z) $ G (4, a) , (48) 
o 

where D^ x \rT, r^, a; b, z) is the amplitude of production of a colorless cc pair which reaches 
a separation f T outside the nucleus. It is produced at the point (6, z) by a color-octet cc 
with separation r' T . The amplitude consists of two terms, 

(f T , 4, a; 6, z) = D {X) (f T , 4, a; b, z) + D [x) (f T , r' T , a; 6, z) . 



Here the first term reads, 

D {X) (f T , 4, a; b, z) = Gg (f T , z + ; r' Tl z) e W ■ d' e^ z 



(49) 



(50) 



where G^i^T-, z+] ?ti z ) is the color-singlet Green function describing evolution of a cc wave 
packet with initial separation r' T at the point z up to the final separation r*r at z + ^ oo. 
This term is illustrated in Fig. [3]a. 

There is also a possibility for the projectile gluon to experience diflractive interaction 
with production of color-octet cc with the same quantum numbers of the gluon at the 
point z\. This pair propagates from the point z\ to z as is described by the corresponding 
color-octet Green function G&) and produces the final colorless pair which propagation is 
described by the color-singlet Green function, as is illustrated in Fig. ||b. The corresponding 
second term in (H91) reads, 



D i 2 x \f Tl f^a;b,z) 



1 

2 



dz x d 2 r'r (f T , z + ; r£, z) 



^ ■d"Gi\f^z-rf r ,z 1 )e lc '^a 3 (rf r ,a)p A (b,z 1 ) , 



(51) 



The singlet, , and octet, GfJ, Green functions describe the propagation of color- 
singlet and octet cc, respectively, in the nuclear medium. They satisfy the Schrodinger 
equations, 



d 



% — G^{f T ,r' T \z,z') 
dz 



2 E G a (1 - a) 



+ V (k \f T ,a) 



(52) 
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a 



Figure 3: The incident gluon can either produce the colorless cc pair with 
quantum numbers of x at the point z (a), or it produces diffractively a 
color-octet cc with the quantum numbers of the gluon at the point Z\ 
which is then converted into a color singlet state at z (b). Propagation 
of a color-singlet or octet cc is described by the Green functions G^ and 
Gfc, respectively. 



with k = 1, 8 and boundary conditions 

Gic } {r T ,f^ z, z' 



8(f T - fy) 



(53) 



The imaginary part of the LC potential is responsible for the attenuation in nuclear 
matter, 

lmV^(f T ,a) = -l-^ k \r T ,a)p A (b,z) , (54) 



where 



a (1) (r T ,a) = crqq(r T ) , 

a (8) (r T ,«) = a 3 (r T ,a) . (55) 

The real part of the LC potential V( k \rj>, a) describes the interaction inside the cc 
system. For the singlet state HeV^(rr, a) should be chosen to reproduce the charmonium 
mass spectrum. With a realistic potential (e.g. see in ||22|| ) one can solve Eq. (^) only 
numerically. Since this paper is focused on the principal problems of understanding of the 
dynamics of nuclear shadowing in charmonium production, we chose the oscillator form of 
the potential [|J, 

ReV^(f T ,a) = n / (a l r2T , , (56) 



2E G a(l 



a) 



where 



a (a) 



2 \J a(l — a) \ii>j , 

0.3 GeV . 



(57) 



m, 
~2 
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The LC potential 



corresponds to a choice of a potential, 

1 



U(f) 



-Q 



in the nonrelativistic Schrddinger equation, 



*(r) = E^(r) , 



(58) 



(59) 



which should describe the bound states of a colorless cc system. Of course this is an approx- 
imation which we are enforced to do in order to solve the evolution equation analytically. 
To describe color-octet cc pairs we should fix the corresponding potential at 



ReV {s \f T ,a) = 

in order to reproduce the gluon wave function Eq. (|7J). 



(60) 



3.4 Numerical estimates 

In order to keep calculations simple we use the approximation Eq. (|i~9"D for the dipole cross 
section which is reasonable for small-size heavy quark systems. Then, taking into account 
Eqs. (|54|) - (|56D we arrive at the final expressions, 



V {k \r T ,a) = ^K^r. 



2 

T 5 



-(1) 



a 4 (a) 



a(l — a) E, 



a 



iC(s) p A , 



K 



(8) 



-iC(s) p A 



9 r 



a 2 + (l-a) 2 



(61) 

(62) 
(63) 



Making use of this approximation and assuming a constant nuclear density pA{b, z) = pa 
the Green functions can be obtained in an analytical form, 



(r T , r' T \ z 2 , Z\) 



b (k) 



2tt sinh(fi( fc ) Az 
f 

x exp < 



2 Tt ■ T, j 

tanh(fi( fc ) Az) ~ sinh(fi( fc ) A; 



(64) 



where 



6 <*> 
Az 



/t( fc ) Eg a(l — a) 
b (k) 



E G a (I — a) 
z 2 - Zx ■ 
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With the oscillator potential we have chosen the wave function of \ has a simple form 

$ x (r T , a) oc (e+ • r T ) exp (-- a 2 (a) r^j , (65) 
which allows to perform analytical integrations over fy and r' T in Eq. (f48|) for the amplitude, 







Here 



(6, 2, a) e iqLZ - 



.(8) 



d*i p A (b, z) e iqLZl W(b, z, z u a) 



-R(b) 



U {b, z, a) 



Fi(Ai 



&W sinh(0i) +a 2 (a) cosh(0i) 



2 ' 



(66) 
(67) 



where 



W(b, z, Z\, a) 



sinh(0 2 )F 1 (A 2 ) 



+ 



&(8) 6(8) sinh(0 2 ) + 92 cosh(0 2 ) 
F 2 (A 2 ) 



fe( 8 ) sinh(0 2 ) + g 2 cosh(0 2 ) 



3 ' 



(68) 



Fi(Ai) = (l+Ai)e Al £i(Ai) - 1 , 

F 2 (A 2 ) = (A' + 4A 2 + 2)e A2 J B 1 (A 2 ) - A 2 - 3 , 

and E\(x) = dt e~ xt /t is the integral exponential function. Further notations are 

4>i = Azi , mueller 

AZx = yjR\-b 2 -Z, 

<p 2 = Az 2 , 

Az 2 = z — z x , 



(69) 



Ai 

A 2 

9 
N 



mt. 



a 2 (a) sinh(0i) + cosh(0i) 



2 6« 

2 

mi 



a 2 (a) cosh(0i) + feW sinh(0i) 
g sinh(0 2 ) + b^ cosh(0 2 ) 



2 6( 8 ) 

2 

2Ai 



9 cosh(0 2 ) + 6( 8 ) sinh(0 2 ) 



1 a.C(jiw)W 



2vr 



m. 
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We define the nuclear transparency for x production as 

a(GA^ X X) 

TrA{x) = Aa(GN^ X X) ■ (70) 

It depends only on the x or projectile gluon energy. We plot our predictions for lead 
in Fig. |j. Transparency rises at low energy since the formation length increases and the 



X 
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Figure 4: Nuclear transparency for x production off lead as function 
of energy of the charmonium, or x 2 (the upper scale). The solid curve 
includes both effects of coherence and formation, while the dashed curve 
corresponds to l c = 0. Since transparency scales in x 2 according to ([71]), 
values of x<i are shown on the top axis. 



effective absorption cross section becomes smaller. This behavior, assuming l c = 0, is 
shown by dashed curve. However, at higher energies the coherence length is switched on 
and shadowing adds to absorption in accordance to Eq. ([44]). As a result, transparency 
decreases, as is shown by the solid curve. On top of that, the energy dependence of the 
dipole cross section makes those both curves for Tr a (E x ) fall even faster. 

Apparently, the nuclear transparency depends only on the x energy, rather than the 
incident energy or x\. It is interesting that it lead to the x 2 scaling. Indeed, the % energy 

Ml 

= — , (71) 

x 2m N x 2 



20 



depends only on x 2 . We show the x 2 scale in Fig. || (top) along with energy dependence. 

We also compare in Fig. |5| the contribution of quark shadowing and absorption (thin solid 
curve) with the nuclear suppression observed at 800 GeV 0. Since data are for Wj Be ratio, 




-0.1 0.1 0.3 0.5 0.7 0.9 

OC p 

Figure 5: Tungsten to beryllium cross section ratio as function of Feyn- 
man xp for J/^ production at proton energy 800 GeV. The thin solid 
curve represents contribution of initial state quark shadowing and finals 
state cc attenuation for x production. The dotted curve includes also 
gluon shadowing. The dashed curve is corrected for gluon enhancement 
at large x 2 (small xf) using the prescription from H38j . The hnal solid 
curve is also corrected for energy loss and for x —* J / ^7 decay. Experi- 
mental points are from the E866 experiment ||/. 

and our constant density approximation should not be applied to beryllium, we assume for 
simplicity that all pA cross sections including pN obey the A a( * XF \ We see that the calculated 
contribution has quite a different shape from what is suggested by the data. It also leaves 
plenty of room for complementary mechanisms of suppression at large Xf- 
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4 Gluon shadowing 



Previously we considered only the lowest |cc) fluctuation of the gluon what is apparently 
an approximation. The higher Fock components containing gluons should be also included. 
In fact they are already incorporated in the phenomenological dipole cross section we use, 
and give rise to the energy dependence of cr^L. However, they are still excluded from nuclear 
effects. Indeed, although we eikonalize the energy dependent dipole cross section the higher 
Fock components do not participate in that procedure, but they have to be eikonalized as 
well. This corrections, as is demonstrated below, correspond to suppression of gluon density 
in nuclei at small x. 

The gluon density at small x in nuclei is known to be shadowed, i.e. reduced compared 
to a free nucleons. The partonic interpretation of this phenomenon looks very different 
dependent on reference frame. In the infinite momentum frame, as was first suggested by 
Kancheli the partonic clouds of nucleons are squeezed by the Lorentz transformation 



less at small than at large x. Therefore, while these clouds are well separated in longitudinal 
direction at large x, they overlap and can fuse at small x, resulting in a diminished parton 
density f3|,g(|. 

Different observables can probe this effect. Nuclear shadowing of the DIS inclusive 
cross section or Drell-Yan process demonstrate a reduction of the sea quark density at 
small x. Charmonium or open charm production is usually considered as a probe for gluon 
distribution. 

Although observables are Lorentz invariant, partonic interpretation is not, and the mech- 
anism of shadowing looks quite different in the rest frame of the nucleus where it should 
be treated as Gribov's inelastic shadowing. This approach seems to go better along with 
our intuition, besides, the interference or coherence length effects governing shadowing are 
under a better control. One can even calculate shadowing in this reference frame in a pa- 



rameter free way (see in |)41| , [28| , ||) employing the well developed phenomenology of color 
dipole representation suggested in P0| . On the other hand, within the parton model one 
can only calculate the Q 2 evolution of shadowing which is quite a weak effect. The main 
contribution to shadowing originates from the fitted to data input. 

In the color dipole representation nuclear shadowing can be calculated via simple eikon- 
alization of the elastic amplitude for each Fock component of the projectile light-cone wave 
function which are the eigenstates of interaction JftJ. Different Fock components represent 



shadowing of different species of partons. The \qq) component in DIS or \q^*) in Drell- 
Yan reaction should be used to calculate shadowing of sea quarks. The same components 
including also one or more gluons lead to gluon shadowing [^. [| . 

In the color dipole approach one can explicitly see deviations from QCD factorization, i.e. 
dependence of the measured parton distribution on the process measuring it. For example, 
the coherence length and nuclear shadowing in Drell-Yan process vanish at minimal X2 (at 



fixed energy) J24], 25], while the factorization predicts maximal shadowing. Here we present 
even more striking deviation from factorization, namely, gluon shadowing for charmonium 
production turns out to be dramatically enhanced compared to DIS. 
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4.1 LC dipole representation for reaction G N — > xGX 



In the case of charmonium production different Fock components of the projectile gluon, 
\{cc)inG) containing a colorless cc pair and n gluons (n=0,l...) build up the cross section 
of charmonium production which steeply rises with energy (see in |[22||). The cross section 
is expected to factorize in impact parameter representation in analogy to the DIS and Drell 
Yan reaction. This representation has the essential advantage in that nuclear effects can 
be easily calculated 
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Feynman diagrams corresponding x production associated 
with gluon radiation are depicted in Fig. [12] of the Appendix D. We treat the interaction 
of heavy quarks perturbatively in the lowest order approximation, while the interaction 
with the nucleon is soft and expressed in terms of the gluon distribution. The calculations 



performed in [Appendix D| are substantially simplified if the radiated gluon takes a vanishing 



fraction a 3 of the total light-cone momentum and the heavy quarkonium can be treated as 
a nonrelativistic system. In this case the amplitude of x G production has the simple form 
Eq. ( p.40| ) which corresponds to the "Drell- Yan" mechanism of x production illustrated in 
Fig. H Correspondingly, the cross section of x production has the familiar factorized form 




+ 




Figure 6: The dominant Feynman diagrams contributing to x production 



(see \Appcndix L\). 



similar to the Drell- Yan reaction p6| , [27| , 
da(GN ^ X GX) 



a 3 



d a 3 



J d 2 s |^ Gx (s,a 3 )| 2 a GG [(1 - a 3 )s,x 2 /a 3 ] , (72) 



where (Jgg{t, x ) — 9/4 aq q (r, x) [see Eq. fll8|) l is the cross section of interaction of a G G dipole 
with a nucleon. \l/(s, a 3 ) is the effective distribution amplitude for the x ~ G fluctuation of 
a gluon, which is the analog to the 7* q fluctuation of a quark, 



^ Gx (s,a 3 ) = Yl / d 2 rdot$f (r, a) $^(f, a) 



x 



r a 3 \ I _ r 

s+-,- I.-- 



a 3 



a 



(73) 



The notations of [Appendix D| are used here, $^ M (r*, a) and 3>cf(r, a) are the qq LC wave 
functions of the x an d gluon, respectively, which depend on transverse separation r and 
relative sharing a by the qq of the total LC momentum. $ x g(s, a) is the LC wave function 
of a quark-gluon Fock component of a quark given by Eq. ( p.36| ). 
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4.2 Gluon shadowing for % production off nuclei 

The gluon density in nuclei is known to be modified, shadowed at small Bjorken x. Corre- 
spondingly, production of x treated as gluon-gluon fusion must be additionally suppressed. 
In the infinite momentum frame of the nucleus gluon shadowing should be treated as a 
GG — > G fusion which diminishes the amount of gluons p^, M. In the rest frame of the 



nucleus it looks very different as Gribov's inelastic shadowing |43[ related to diffractive gluon 
radiation. This frame seems to be more convenient way to calculate gluon shadowing, since 
techniques are better developed, and we use it in what follows. 

The process G N — > x X considered in the previous section includes by default radiation 
of any number of gluons which give rise to the energy dependence of the dipole cross section 
in© 

Extending the analogy between the reactions of x G production by an incident gluon and 
heavy photon radiation by a quark to the case of nuclear target one can write an expression 
for the cross section of reaction G A — > x G X in two limiting cases: 

(i) the production occurs nearly instantaneously over a longitudinal distance which is 
much shorter than the mean free path of the xG pair in nuclear matter. In this case the 
cross sections on a nuclear and nucleon targets differ by a factor A independently of the 
dynamics of x G production. 

(ii) The lifetime of the x G fluctuation, 



tr 



2 E G 

Kg 



(74) 



substantially exceeds the nuclear size. It is straightforward to replace the dipole cross section 
on a nucleon by a nuclear one ^B], then Eq. ( [72] ) is modified to 



da(GA -> X GX) 
d (In a 3 ) 



2 / d 2 bd 2 s Itf 



1 — exp 



• (75) 



In order to single out the net gluon shadowing we exclude here the size of the cc pair 
assuming that the cross section responsible for shadowing depends only on the transverse 
separation s. 

A general solution valid for any value of t c is more complicated and must interpolate 
between the above limiting situations. 

(iii) In this case one can use the methods of the Landau-Pomeranchuk-Migdal (LPM) theory 
for photon bremsstrahlung in a medium generalized for targets of finite thickness in |37], f2"8fl . 
The general expressions for the cross section which reproduces the limiting cases t c — > (i) 
and t c — > oo (ii) reads, 
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where 

£(22,21) 



d s x d s 2 *g x (^2, "3) o- GG (s 2 , x 2 /a 3 ) G(s 2 , z 2 ; si, 2 X ) a GG {s x , x 2 /a 3 ) *g x (^i ? "3) 

(77) 

Here the Green function G(s*2, 22; si, z\) describes propagation of the \G P & i r which starts 
with transverse separation si at the points with longitudinal coordinates z\ and ands up at 
z 2 having separation si. It satisfies the Schrodinger type equations, 



dG(s 2 ,z 2 ; si,zi) 
dz 2 

dG(s 2 ,z 2 ; si,2i) 



A 



« 12 



In the limit z 2 - 



2E G a 3 {\ 
A 



"3 



81 



-ia GG (s 2 ) p A (B,z 2 ) 
-ia GG (si) p A (B,Zi) 



z\ 



2E G a 3 {\ - a 3 ) 
it should satisfy the condition 
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G(s 2 ,z 2 ; si, zi) ; (78) 
G(s 2 ,2 2 ;si,2i) . (79) 

(80) 



A full calculation of the cross section of associated \ G production off nuclei Eq. flTE|) 
employing the solutions of equations (|77|) - (p0| ) with realistic shapes of the dipole cross 
section and nuclear density can be done only numerically, and is still a challenge for com- 
puting. However, the problem can be essentially simplified if the following approximations 
are done, 



PA (b,z) = q(r a -VvTz^) ; 

vgg(s,x) w C g (x) s 2 . 



(81) 



The former is rather realistic for heavy nuclei, while the latter needs special care to be 
adjusted to realistic calculations, since the value of factor C G (x) correlates with the mean 
separation (s) which depends on the process. 

Calculations of the nuclear cross section with the realistic dipole cross section which 
levels off at large separations is complicated only for the transition region of l c ~ R A . In 
the limit of l c — > 00 one can easily perform calculations with any form of the dipole cross 
section and adjust the C G (x) to reproduce the cross section of reaction G A — ► \ G X. Such 
a prescription guarantees the correct endpoint behavior at t c ^> R A and t c — > 0, then we 
expect the transition region should not be very wrong either. 

To single out the correction for gluon shadowing one should compare the cross section 
Eq. ([7E|) with the impulse approximation term in which absorption is suppressed, 
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For further calculations and many other applications one needs to know gluon shadowing 
as function of impact parameter which is calculated as follows, 
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where da(GA — > \GX) / (Pbdotz is the function integrated over <i 2 5 in (|76|) . The results of 
calculations for the 6-dependence of gluon shadowing (|83| ) are depicted in Fig. |7| for different 
values of x<i as function of thickness of nuclear matter, L - 



R\ — b 2 . The results confirm 
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Figure 7: Gluon suppression as function of thickness of nuclear matter 
with constant density pA = 0.16 fm~ 3 . 



the obvious expectation that shadowing increases for smaller %2 and for longer path in 
nuclear matter. One can see that for given thickness shadowing tends to saturate down to 
small X2, what might be a result of one gluon approximation. Higher Fock components with 
larger number of gluons are switched on at very small Xi- At the same time, shadowing 
saturates at large lengths what one should have also expected as a manifestation of gluon 
saturation. Note that at large x 2 = 0.03 shadowing is even getting weaker at longer L. 
This is easy to understand, in the case of weak shadowing one can drop off the multiple 
scattering terms higher than two-fold one. Then the shadowing correction is controlled 
by the longitudinal formfactor of the nucleus which decreases with L (it is obvious for 
the Gaussian shape of the nuclear density, but is also true for the realistic Woods-Saxon 
distribution). 

The next problem we face is how to correct our previous calculations for the calculated 
gluon shadowing. The standard parton model prescription is to multiply the cross section 
of charmonium production by Rq{x2). This may be correct only if the process is so hard 
that no nuclear effects except gluon shadowing exists. For instance, this is the case in DIS 
for highly virtual longitudinally polarized photons. However, as we calculated and Fig. £| 



26 



demonstrates, multiple interaction of the cc pair cannot be ignored, and the parton model 
prescription is invalid. To get a correct result one should perform a full calculation of nuclear 
effects which involve all Fock components. This is still a challenge, meanwhile one should 
look for a reasonable approximation. 

Our approximation of the lowest Fock state containing only one gluon corresponds in 
term of Regge theory to Pomeron fusion nW — > TP (n=2,3...), while transitions nTP — > mTP 
(m > 2) are missed. This is explained pictorially in Fig. |[ Although inclusion of the 




Figure 8: Multiple rescattering of a gluonic component of a hadron in 
the target rest frame. In t-channel it is interpreted as multi-Pomeron 
fusion n IP — > IP. 

multi-gluon Fock states is problematic, there is a simple and intuitive way to include an 
essential part of missed corrections. Intuitively one can say that the projectile high energy 
partons experience multiple interactions not with target nucleons, but with gluons at small 
X2- Eikonal (Glauber) approximation assumes that the number of gluons is proportional 
to the number of nucleons. This is why the eikonal exponent contains a product of the 
interaction cross section and nuclear thickness, <jT a (B). It is not true any more if gluon 
fusion is at work. The gluon density decreases (the more, the larger the T A is) and one can 
easily take it into account renormalizing the nuclear density, 

p A {b, z) p A (b, z) = R G (x 2 , b) p A {b, z) . (84) 

In this way we do account for the missed nTP — > mTP transitions. Indeed, m-fold interaction 
associated with mTP exchange is not are not in sequence, this would lead to the famous 
AFS (Amati-Fubini-Stangelini) cancelation. The m Pomerons correspond to simultaneous 
development of m parton ladders. In our approximation each of this Pomeron ladder is 
suppressed by fusion precess which is exactly n'TP — > IP in this case. By eikonalizing [see 
e.g.Eq. (|44|)1 the result of fusion we correctly involve the higher orders of n'IP — > IIP, as is 
illustrated in Fig. ^. 

Thus, we correct our results for for gluon shadowing renormalizing the nuclear density 
by means of (jS4|). It must be done in the imaginary part of the light-cone potential (|54"D in 
the evolution equation Eq. ( |52"D and also in Eqs. (|^) and fl5T~|). The numerical results are 
depicted by dotted curve in Fig. |5[ Apparently, gluon shadowing is stronger at small x 2 , i.e. 
at large xf- 
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Figure 9: Multiple rescattering of a projectile hadron (e.g. cc) via multi- 
Pomeron interactions diminished by nlP — > IP fusion. 



4.3 Antishadowing of gluons 

Nuclear modification of the gluon distribution is poorly known. There is still no experimental 
evidence for that. Nevertheless, the expectation of gluon shadowing at small x is very solid, 
and only its amount might be disputable. At the same time, some indications exist that 
gluons may be enhanced in nuclei at medium small £2 ~ 0.1. The magnitude of gluon 



antishadowing has been estimated in [44J] assuming that the total fraction of momentum 
carried by gluons is the same in nuclei and free nucleons (there is an experimental support 
for it). Such a momentum conservation sum rule leads to a qluon enhancement at medium 
x, since gluons are suppressed in nuclei at small x. The effect, up to ~ 20% antishadowing 
in heavy nuclei at x ~ 0.1, found in |44| is rather large, but it is a result of very strong 
shadowing which we believe has been grossly overestimated (see discussion in ||). 

Fit to DIS data based on evolution equations performed in |3S| also provided an evidence 
for rather strong antishadowing effect at x ~ 0.1. However, the fit employed an ad hoc 
assumption that gluons are shadowed at the low scale Qq exactly as ^(x, Q 2 ) what might 
be true only by accident. Besides, in the x distribution of antishadowing was shaped ad hoc 
too. 

A similar magnitude of antishadowing has been found in the analysis |45j of data on 



Q 2 dependence of nuclear to nucleon ratio of the structure functions, F^(x, Q 2 )/ F^ (x, Q 2 ). 
However it was based on the leading order QCD approximation which is not well justified 
at these values of Q 2 . 

Although neither of these results seem to be reliable, similarity of the scale of the pre- 
dicted effect looks convincing, and we included the antishadowing of gluons in our calcu- 
lations. We use the shape of Xi dependence and magnitude of gluon enhancement from 
[38] . 



5 Energy loss corrections 

The mechanism of J/\l/ suppression at large xp due to initial state energy loss of the pro- 
jectile partons in nuclei was first suggested back in 1984 ||23|| . As soon as the incoming 
hadron interacts inelastically in the nucleus, the projectile partons (only those which are 
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sufficiently soft to be resolved by the soft interaction) lose coherence and start losing energy 
for hadronization. If the coherence time of charmonium production is shorter than the mean 
internucleon separation in a nucleus, as it was assumed in |p3| , the projectile parton (gluon) 
responsible for charmonium production arrives at the production point with a diminished 
energy Eq — AEq- In order to produce a charmonium with the same energy (fixed by the 
measurement) as on a proton target the projectile gluon must have an excess in the initial 
momentum which leads to a suppression 
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g h (xi + Axi 
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where g (x) is the gluon density in the incoming hadron h and 
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However, with some maybe small probability Wo, the incoming hadron may have no 
initial state inelastic interactions and produce charmonium in the very first interaction 
without any energy loss. Thus the nuclear suppression acquires two contributions 
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Here the energy loss AEq depends on the longitudinal distance z 2 — z\ between the first 
inelastic interaction and the production point of the charmonium. The probability of no 
initial state interaction reads, 
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The second term in ( |87| ) is normalized to the probability W\ 
case of zero energy loss [Rq = 1), 



1 — Wo of interaction in, the 
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(89) 



In order to calculate nuclear suppression (j87D we need to know the function dEo/dz 
which the rate of energy loss of a gluon propagating through nuclear matter. Energy loss 
initiated by the inelastic collision continues with a constant rate like it were in vacuum. One 
comes to this conclusion both in the color-string model p3| , [46[ or treating hadronization as 
gluon bremsstrahlung 
have a rate 



|48fl . This so called vacuum energy loss was supposed in p3| to 
-dEc/dz = 3GeV/ fm leading to a rather good description of data from the 
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NA3 experiment JT|. A substantially larger value —dEc/dz = 5 GeV/ fm was advocated in 



[ 35| for the intermediate color-octet state in photoproduction of J/^f off nuclei. 

In fact, there is no controversy here, because in the case of hadroproduction of char- 
monium it is more appropriate to rely upon the energy loss of a quark, rather than a 
gluon, propagating through a medium. Indeed, it was found in || that light-cone gluons 
are located at small transverse separation ~ 0.3 fm from the valence quarks. This is dic- 
tated by data for large mass diffraction which corresponds to diffractive gluon radiation and 
was interpreted in as a result of a strong nonperturbative interaction of the light-cone 
gluons. This observation is in a good accord with other models for the nonperturbative 



QCD effects (see discussion in [^9j). Thus, only a semi-hard interaction can resolve a gluon 



in the constituent quark, while a soft interaction responsible for inelastic interactions of 
the incident proton in the nucleus do not see the gluon. The whole constituent quark, 
rather than the gluon, is hadronizing and loosing energy propagating through the nucleus. 
Therefore, we should expect a rate of energy loss similar to what is observed for Drell-Yan 



lepton pair production off nuclei. The recent analysis p3| of data [12| for nuclear sup- 
pression in Drell-Yan reaction revealed for the first time a nonzero energy loss for quarks, 
—dEg/dz = 2.32 ± 0.52 ± 0.5 GeV/ fm. This is compatible with the value 3 GeV/ fm used 
in @, 0. 

Thus, one can calculate the nuclear suppression of charmonium production caused by 



energy loss in the same way as for Drell-Yan reaction [24, 25 



oo 1 

fdLW x {L) J dx q F q h (xq) daf{x\)/dx\ 

Rf/N^i) = w + 1 (Xq)mm . (90) 

fdx q F^Xq)daf{x\)/dx\ 

xi 

Here F^{x q ) is the quark distribution function in the incident hadron, and x q and x\ = 
Xi/Xg are the fraction of the light-cone momentum of the incoming hadron h carried by the 
quark and the fraction of the quark momentum carried by the cc, respectively. The lower 
integration limit is given by {x q ) min — x\ + AE/E h , and x\ = x\/ [x q — AE/E h ). The cross 
section cr^ixf) of \ production by a constituent quark, qN — > \ X is assumed to behave 
as oc (1 — xf) 2 ' 5 in order to reproduce the observed 1 — xi) 5 distribution of J/\l/ production. 

At first glance the case of long coherence length is more complicated since the cc fluc- 
tuation is produced long in advance and propagates through the nucleus rather than the 
projectile gluon. However, the soft interaction responsible for the first inelastic h — N colli- 
sion (see above) does not discriminate between the gluon and the color-octet cc. The same 
is true if instead of the light-cone representation one uses the equivalent description of co- 
herence via diffractive transition G —> cc with no color exchange. Therefore, this inelastic 
interaction of the projectile hadron and energy loss modify the ratio (^) in the same way 
as is described above for the limit of short coherence time. Thus, Eq. (^) is valid in general 
case for any values of coherent and formation lengths. 

Further, multiple interactions of a parton in a nuclear medium leads to a broadening 
of the parton's transverse momentum and enhanced gluon bremsstrahlung. The rate of 
associated induced energy loss rises linearly with the length of the path in nuclear matter 



| 50| , but is a rather small correction to the vacuum energy loss. 
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The induced energy loss fluctuates and it becomes an important effect towards the 
kinematical limit X\ — > 1 where no gluon can be emitted. This of course suppresses the 
production rate of charmonium, and more on a nuclear than on a proton target. Indeed, 
the projectile gluon experiences broadening of the transverse momentum due to multiple 



interactions (see in [51] connection between px broadening and induced energy loss). As a 
result gluon radiation becomes more intensive. The relation between induced energy loss and 
nuclear broadening of the mean transverse momentum squared of the quark was established 
in [f)0|| , 

AE = ^A(p 2 T )L, (91) 

where the broadening of (pj) is proportional to the length L of the path in nuclear matter. 
The p\ broadening for J/\l/ was measured in the E772 experiment |Q to be rather small 



A(pt) ~ 0.5 GeV 2 even for tungsten. Therefore the effective rate of induced energy loss 
in fl9~T| ) is about 0.37 GeV/ fm, nearly an order of magnitude smaller than the the value 
of vacuum energy loss we use. We assume that the induced energy loss is incorporated 
effectively in dE/dz = —3 GeV/ fm we use in our calculations. 



6 Modification of the cc^-distribution by % — > J/\I>7 
decay 

No data is still available for x^-distribution of produced x? but only for J/ty. About 40% 
of them originate from the x ~^ decay and have momenta smaller than that of chi. 

Therefore we should correct the x f distribution calculated in previous section to compare 
with data for indirect J/^> . It would be incorrect, however, to assume isotropic decay of x 
since it is produced polarized. 

The structure of the amplitude of decay of the tensor meson (x 2 ) to two vector mesons 
with masses m\ and m 2 is fixed by the gauge invariance, 



A(T -> V x + V 2 ) oc hT 
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where e\ )2 are the polarization vectors of Vi j2 , 

eiki = e 2 k 2 



; 



(93) 



(94) 



p is the 4-momentum of the X2', ^1,2 are the 4-momenta of V\ and V 2 respectively; h^ v = h Ufl 
is the X2 polarization tensor satisfying the conditions, 



p^ 



; 
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(95) 
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In the rest frame of the \2 the components of the polarization tensor h 00 = h° l = h l ° = 0. 
Other components for the state with the projection m = 2 are, 

h ik {m = 2) = t\e% , (96) 

where 

? + = -^r- < 97 > 

The angular distribution of the decay products relative to the momentum of Xi nas the 
form, 

oc (1 + /3x cos 2 6)(l + f3 2 cos 2 6) , (98) 

a cosy 

where [3^ are the velocities of the produced Vi }2 - In the case of V\ — J/^, V 2 — 7 we have 
/3i « 1/9', ft = 1. 

7 Comparison with available data and predictions for 
higher energies 

First of all, we checked with data at the energies of the SPS where the coherence length is 
rather short and the main effects are absorption and energy loss. With (|9~0"D we calculated the 
xf dependence of the exponent a(xp) describing the A-dependence of the J/^/ production 
rate, oc A a ( XF ^ at 200 GeV. We also corrected it for gluon enhancement and x ~^ J/^"f 



decay. The results are compared with data in Fig. [II]. The observed steep fall of a is 
well reproduced, while at small xp our calculations seem to underestimate the data from 
the NA3 experiment. In fact, our calculations for energy loss effect are not trustable at 
small xf- Indeed, the suppression caused by the shift of energy depends on how steep 
the x\ dependence of the cross section qN — > %X is, the steeper it is, the stronger is the 
suppression. The chosen parametrization da^ N j d x\ oc (1 — x\) 2 ' 5 is valid only at large 
x\ — > 1. This DY cross section in iViV collision must be flat around xp = 0, what makes 
the suppression by energy loss weaker at small xp. 

At higher energy, 800 GeV the dynamics of J/\l/ suppression is more complicated and 
includes more effects. Now we can apply more corrections to the dotted curve in Fig. |5] 
which involves only quark and gluon shadowing. Namely, inclusion of the energy loss effect 
and decay x —> J/^ 1 leads to a stronger suppression depicted by dashed curve. Eventually 
we correct this curve for gluon enhancement at x 2 ~ 0.1 (small xp) and arrive at the final 
result shown by thick solid curve. 

Since our calculation contains no free parameters we think that the results agree with the 
data amazingly well. Some difference in the shape of the maximum observed and calculated 



at small xp may be a result of the used parameterization 0] for gluon antishadowing. 
We think that it gives only the scale of the effect, but neither the ad hoc shape, nor the 
magnitude should be taken literally. Besides, our calculations are relevant only for those 
J/^s which originate from x decays which feed only about 40% of the observed ones. 

At higher energies of the RHIC and LHC the effect of energy loss is completely gone and 
nuclear suppression must expose the x 2 scaling. Much smaller x 2 can be reached at higher 
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Figure 10: xf dependence of nuclear suppression of J/^fr production at 
200 GeV in terms of the exponent a characterizing the A a dependence 
of J/ty production. Energy loss Eq. (jffij), absorption, quark shadowing 
and gluon enhancement are included in the calculations. The round data 
points are from |2|/, the square point is from J53| / . 



energies. Our predictions for proton-gold to proton-proton ratio is depicted in Fig. [TT| One 
can see that at xp > 0.1 shadowing suppresses charmonium production by nearly an order 
of magnitude. 

We can also estimate the effect of nuclear suppression in heavy ion collisions assuming 
factorization, 

Rab(xf) = R p b{x f ) R- p a(-xf) ■ (99) 

Our predictions for gold-gold collisions at ^fs = 200 GeV are shown by the bottom curve in 
Fig. |ll[ Since factorization is violated this prediction should be verified. 

8 Conclusions and outlook 

The long standing challenge of explanation of the steep xf dependence of charmonium sup- 
pression observed in pA collisions becomes an appealing problem with advent of RHIC and 
LHC colliders. A good understanding of the underlying mechanisms is especially important 
since charmonium suppression is one of the major probes for the quark-gluon plasma forma- 
tion in heavy ion collisions. The conventional mechanisms of charmonium suppression used 
as a baseline for new physics in experiments at SPS CERN must be essentially revised. This 
demands development of new theoretical tools, and first of all an explanation for already 
collected data for charmonium suppression. 

Coherent effects or shadowing become the dominant effect which governs charmonium 
suppression at RHIC. So far no theoretical tool has been available to deal with this phe- 
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Figure 11: Nuclear suppression of J/ty production in proton-gold col- 
lisions at y/s = 200 GeV as function of xp (the upper curve) and in 
gold-gold collisions (bottom curve). Effects of quark and gluon shadow- 
ing and gluon antishadowing are included. 



nomenon. The main result of present paper is the development of the light-cone dipole 
approach describing coherence effects at any energy. It can be structured as follows. 

• Final state absorption of the produced colorless cc in nuclear matter is accomplished by 
initial state shadowing at high energies when the coherence time Eq. (0) is sufficiently 
long. A light-cone Green function formalism is developed which describes propagation 
of fluctuating color octet (singlet) cc wave packets through nuclear matter. This 
approach includes both the effects of formation and coherence. It interpolates between 
the low energy limit, tf <C Ra, where only final state absorption causes charmonium 
suppression, and the high energy limit, t c ^> Ra where the color octet cc fluctuation 
is developed long in advance and propagates through the whole nuclear thickness and 
is converted to a color singlet. This is an analog to shadowing of c quarks in DIS. Its 
existence is a deviation from QCD factorization which is quite strong, an about 50% 
effect on heavy nuclei. 

• At high energies fluctuations are overwhelmed by gluons which are the source of rising 
energy dependence of all cross sections. Such fluctuations are rather heavy and their 
lifetime is much shorter than that of a cc pair. Only at very high energies these 
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gluonic fluctuations live sufficiently long to participate in multiple interactions which 
lead to shadowing. This gluon shadowing turns out to be very strong and becomes 
the main source of charmonium suppression at RHIC. Again QCD factorization is 
strongly violated. Namely, gluon shadowing for colorless cc production turns out to 
be much stronger than in DIS of open charm production. 

• Both mechanisms, shadowing of c quarks and gluons depend only on Xi- Although 
available data at medium high energies strikingly contradict £2 scaling, it comes from 
energy loss effects. Indeed, modification of the x\ distribution by the energy loss de- 
pends both on x\ and incident energy. Accidentally, data demonstrate an approximate 
xp scaling in the energy interval 200 — 800 GeV which our calculations confirm. It is 
a result of interplay between energy loss and shadowing: the former effects vanishing 
with energy, while the latter is rising. Exact X2 scaling is predicted at high energies 
(well beyond 1 TeV) of RHIC and LHC We predict a very steep variation of charmo- 
nium suppression around xf = at RHIC. In fact, it can also be obtained directly from 
the E866 data neglecting energy loss corrections at 800 GeV . For gold-gold collisions 
we predict a peculiar narrow peak at xf = 0. 

As a further development we are going to extend this approach to the case of directly 
produced J/^f and We do not expect the results to be very different what what we found 
in this paper. Indeed, the kinematics of decays x ~^ J/^'y and | cc) 8 — > J/ty G is similar up 
to the difference in the invariant masses of the cc pairs. However, we demonstrate in Sect. |5] 
that the decay has practically no effect, and the xp distributions of x an d J/^f are nearly 
identical. 

Another source of direct production is hidden in the nuclear effects for x production. 
Namely, the cc pair traveling through nuclear matter and experiencing multiple interactions 
changes its spin-orbital structure, therefore it also feeds the \& channel. We eliminated this 
channel in Sect. || projecting the cc on the state with the quantum numbers of x- If to 
project the cc on the \I> one gets a cross section comparable with data. Therefore, this 
channel of direct \l/ production is important and should be studied. 

It is still a challenge how to apply the developed approach to nucleus-nucleus collisions. 
The relation Eq. fl99|) between nuclear effects in pA and AA collisions is just an approxima- 
tion based on QCD factorization which as we saw is badly violated. We also hope to make 
a progress in this direction. 

Note added: Soon after this paper has been submitted to the Los Alamos bulletin 



board, another paper was put on the server |57j which focuses at the same problem of 
charmonium production at the energies of RHIC. Nuclear suppression at RHIC is predicted 
only at very small xf < 0.084 assuming X2 scaling and employing the E866 data ||. At 
larger xp nuclear effects are calculated within model fitted to the E866 data. The 
model has enough free parameters to be successful with available data, but it completely 
misses the coherence effects, in particular gluon shadowing, which are the main source of 
J/ty suppression at high energies. Therefore, extrapolation to higher energies beyond the 
X2 range of the E866 data should not be trusted. 

On the other hand, we do not rely on either X2 scaling (which is still not at work at 
800 GeV, see Fig. |), or the E866 data, but directly calculate nuclear effects at the energy of 
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RHIC. We compare our parameter free calculations with the E866 data only to demonstrate 
agreement. 
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Appendix A LC wave functions of the % states 

We are going to establish relation between the wave functions of the 3 Pj states (xj mesons) 
in the LC and rest frames. First, the LC wave function in coordinate space, 
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d 2 p T e~ l ^ T ^f (p T ,a) 



is normalized as, 



(A.l) 



1 -L 

52 J da J d 2 r T |$f (r T ,a)| 2 = £ j da J d 2 p T |$f (p T ,a 



MM o 



MM o 



(A.2) 



The problem of Lorentz boosting of a quarkonium wave function from the rest frame, where 
it can be found as a solution of the Schrodinger equation, to the LC frame still has no well 
grounded solution. A nonrelativistic two-body system cc in its rest frame corresponds to the 
full series of Fock states, which include gluons and the sea, in the LC frame. To construct 
the wave function for the lowest Fock component cc we use the popular recipe suggested in 
PJ (see also in pR G2|), 



' dp L (p T , a) 
da 



P=Pt+Pl™ ' 



(A.3) 



where 



p L = (a -1/2) M{p T ,a) 



M 2 (p T ,a) 



m 2 + p\ 
ail — a) 



The rest frame wave function ^^(p) for a state with the total angular momentum J and 
projection J z on the direction n is a combination of different spin and orbital momenta, 



(A.4) 



where the orbital part ipi (p) is a Fourier transform of the spatial part of the wave function, 



(2vr) 3 / 2 



rf 3 re^> ljm2 (f) . 



(A.5) 
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The latter, as usual, can be represented as a product of the angular and radial r-dependent 
parts, 

ii>i, ma (r f )=Yi,m(f f /r)R 1 (r) , (A.6) 

It is normalized as, 

/ d 3 r |^, m2 (f)| 2 = J drr 2 R\{r) = 1 . (A.7) 
The orbital part has the form, 



Y i,m 2 (r/r) = 

where the polarization vectors are defined as, 




3 6rri2 ' ^ 



47T r 



(A.8) 



e+ = -(e x + ie y )/V2 
e = e z = n 

k e+ = (e x -ie y )/y/2 



m = 1 
m = 
m = — 1 



(A.9) 



They satisfy the relations, 



e+ • e_ = -1 , 

— » — » — # — # — # — # — * — * r\ 

e + • e + = e_ ■ e_ = e + • eo = e_ ■ eo = U . 



The spin part reads, 



^mi ~~ ^2 ^ ° ' emi ^ ' 



(A.10) 

(A.ll) 
(A.12) 



where 

and spinors £ in the charmonium rest frame are related to those, 77, in the LC frame by the 
Melosh spin rotation transformation, 



$ = R(-p T ,l-a) V § 



The rotations matrix reads, 



R(pr, ol) 



m c + a M — ia • (n + pr) 



yj \m c + a M) 2 + p 2 T 
Then spin part of the LC wave function takes the form, 



(A.13) 



(A.14) 



where 



A mi (pr, a) = RHpt, «) (<? • e mi ) a y R*(p T , a) a y 1 = B mi + C mi , 
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(A.15) 
(A.16) 



a 



B 



mi 



rriT 



m c a ■ e mi + i(e mi x n) ■ p T 



mi 



(e mi • Pt){o ■ Pt - VlQ ■ n) + (e mi ; n)(p L a ■ p T + p 2 T a ■ n) 

m T (2m c + M) 



and m T = \Jrn 2 + p\. 



(A.17) 
(A.18) 



Appendix B The ol distribution of the charmonium 

wave function 



The a-dependence of the charmonium LC wave function picks at a 
related to the mean longitudinal velocity (vj) of the quarks, 



a — 



<P|> 
Am 2 c 



(vl), 



1/2 and has a width 



(B.1) 



which is a rather small for nonrelativistic systems. 

The relative cc velocity can be estimated applying the Virial theorem to the realistic 
power-law potential V(r) = — Vo+Vi(r/ro) A which describes well the properties of charmonia 
| 56| . Then, the mean kinetic energy (Ek) = m c (v 2 ) is related to the power A, 



(Ek) 



A 



2 + A 



Mm> - 2m c + V 



(B.2) 



where My is the mass of the charmonium. Using A = 0.1 and Vo = 8GeV we get 
(v 2 ) = 0.21. Since Vo is so large the mean velocity turns out to be nearly the same for all 
charmonia {J/^f, X- Numerical calculations [^] confirm this, and also show that the 



results is the same for all realistic potentials. This is not surprising since all those potentials 



look pretty similar (see in 22 



Note that the x states with projections J z = ±1 are not produced since the gluons are 
transversely polarized. The X2 with J z = is not produced at all in the nonrelativistic limit 
v —>■ and is suppressed by the tiny factor v A . Therefore, only the X2 with projection 
J z = ±2 is of practical interest. 

For the S-wave state the mean longitudinal speed squared in Eq. (B.l|) would be (vj) = 
(v 2 )/3. However, for the P-wave state with J z = ±2 we are interested in, the corresponding 
factor is even smaller, (vj) = (v 2 )/5, leading to a very small value of ((a — 1/2) 2 ) = 
(v 2 )/20 = 0.01. 



Appendix C Landau- Yang theorem: manifestation of 

relativistic effects 

Obviously, the spin rotation formalism is rather complicated, and one might hope that 
for nonrelativistic systems like charmonium the corrections are small since R(pr, ot) — > 1 
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when v — > 0. In many cases it is indeed true, for instance the cross section of charmonium 
photoproduction increases by 30% (not a negligible correction either) when the spin-rotation 
effects are included [p2 ]. In some cases, however, the spin rotation is very important. 
Photoproduction of the 25" state for instance, is enhanced by factor 2 — 3 [Q. In the 
case under discussion of hadroproduction of x states the spin rotations also plays a principal 
role. Its suppression leads to a rather intensive production of 3 Pi even by two massless 
gluons forbidden by the Landau- Yang theorem. We demonstrate here that inclusion of the 
spin-rotation resolves this controversy. 

Since production of any state with projection J z = ±1 by two on-mass-shell gluons is 
forbidden we should consider only the xi state with J z = 0. Its wave function in the rest 
frame is a linear combination of the spin and orbital parts with known Clebsch-Gordan 
coefficients, 



V2 

Substituting here Eqs. ( [A.15|) -( |08| ) we get, 

1 



(s_ 1 + c_ 1 )Vi, + i(p) 

This expression can be simplified using relations, 

ipi, m (p) = -i e m ■ p Ri (p) 

where 



Ri(p) 



dr 



r 



J 3/2 (pr) Ri(r) 



(C.2) 



(C.3) 



(C.4) 



and also, 



a 



±i 



2 (e±i ■ p)(a ■ pr -p L a ■ n) 



m T (2m c + M) 



(C.5) 

The terms proportional to C cancel and the wave function of the state with J = 1, J z = 
gets the simple form, 



^ { 



2 nrix 

m c a ■ e _ + z(e_ x n) • pr 



m c a ■ e + + i(e + x n) ■ p T 
e+-p T } Ve ■ 



p T 



(C.6) 



Using this expression we eventually arrive at the LC wave function which has the form, 



$^(r T ,«) 



m c a ■ e + + i(e+ x n 
m r a ■ e_ + z(e_ x n) • V 



e_ ■ tt 



where 



H(r T ) 



rriT tt 



dp T p 2 T Ji(p T r T ) 



dpL 
da 



Ri(p) ■ 



(C.7) 



(C. 



39 



Eq. (|C.7|) can be further simplified taking into account relations, 



(a ■ e + ) e_ — (a • e_) e + = o x (e_ x e+) 
e + x e _ = in , e + • e _ = —1 
e + x n = ie + , e_ x n = ze_ , 

(e+ -f r )(e_ -r T ) = -y . 



The final result reads, 



(C.9) 



$^(f r , a) = -- t# (a x n) ■ r T H(r T ) + 



1 d 



r 2 T H(r T ) 



(CIO) 



Now let us consider the function, 



A = E / da / d V ^f(r T , «)* (e T ■ r T ) $g (f T , 



a 



(C.ll) 



which coincides with the amplitude of Xi production in gluon-nucleon collision up to the 
factor JC(s)/8, if the approximation of dipole cross section aq q (rT, s) = C(s) r T is applied. 



This amplitude can be represented as, 



A = A 1 + A 2 



(C.12) 



where 



A 1 = im 2 c J da J d 2 r T (e x n) ■ r T H{r T ) (e • r T ) K (er T ) ; (C.13) 



Ao 



i da d 2 rT 



1 d 



Tt drT 



r% H(r T ) (e T ■ r T ) (e x n) • V Ko(er T ) 



-ie I da I d 2 r T \ -f— \r T H(r T )} (e T • r T ) (e x n) ■ f T i^i(er T ) , (C.14) 
r T drT 



using the general property of the spinors n, 



E(^ T «€)*(^ f &^)=Tr(atS) , 



(C.15) 



where a and 6 are any two-dimensional matrices. 

Integration over azimuthal angle in Eqs. ( |C.12|) -( C.14 ) can be be performed using rela- 
tion, 

J d<pr Ti r Tk = Tr5 ik r T , (C.16) 
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and we arrive at the following common form for A\ and A 2 , 



Ai t 2 = in (ex ex) ■ n J da Ixp,{oi) 



(C.17) 



where 



h(a) 



h(a) 



J dr T r^K (er T ) H(r T ) , 
d 



o 

oo 



e / drj< rj* K\ (ery) 



drT 



r 2 T H(r T ) 



Integration in the last equation can be performed by parts using the relation, 

d 



drT 



The result reads, 



tt Ki(erx) = —£r T K (er T ) 



J drx K (erT) H(tt) 



(C.18) 
(C.19) 



(C.20) 



(C.21) 



Thus, we conclude that in the limit of on- mass-shell gluon, (e — > m c ), I2 — > —h, arid the 
amplitude (|U.12 ) of xi production vanishes in accordance with the Landau- Yang theorem. 
This is not an obvious result, indeed, the procedure of of Lorentz boost to the infinite 
momentum frame for the quarkonium wave function is ill defined, as was mentioned above. 
It is important that it survives such a rigorous test and recovers the Landau- Yang theorem 
in the light front reference frame. 



Appendix D Gluon radiation process: 

G a + AT — (cc)+ + G b + X 

Nuclear shadowing of gluons is treated by the parton model as glue-glue fusion in the infinite 
momentum frame of the nucleus. On the other hand, in the rest frame it is expressed in 
terms of the Glauber like shadowing for the process of gluon radiation, 

G a + N-^ (cc)f + G b + X , (D.l) 

and production of cc pair with positive C-parity in a color-singlet state. G a ^ are gluons in 
color states a and b. Switching to impact parameter representation one can easily sum up 
all multiple scattering corrections for this reaction on a nucleus, since they take a simple 
eikonal form Besides, one can employ the well developed color dipole phenomenology 
with parameters fixed by data from DIS. 

The amplitude of the process Eq. ( D.l ) is described in Born approximation by the set 



of 12 Feynman graphs depicted in Fig. O. The produced cc pair is assumed to be colorless 
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© @ 



Figure 12: Born graphs contributing to cc pair production accompanied 
with radiation of a gluon. The cc is in a color-singlet and C-even state. 
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and to have positive C-parity. Note that the two sets of diagrams, 1-10 and 11-12, are gauge 
invariant over the t-channel gluon, i.e. the two contribution to the cross section are infra-red 
stable. However, they correspond to different Fock components of the incident gluon in the 
light-cone approach. The first set (1-10) corresponds to a fluctuation containing a colorless 
cc and a gluon (see below), white the second one (11-12) describes a fluctuation with a 
colored cc and a gluon. In the latter case the cc has to interact with the target to become 
colorless, in the same way as for the G — *> cc Fock component considered above in Sect. |[ 
The gluon does not interact at all at the Born level, therefore these graphs 11-12 do not 
contribute to gluon shadowing. 

The amplitude of a colorless cc production corresponding to graphs 1-10 has the following 
structure, 

10 

A^ = Y / A? fl , (D.2) 
i=i 

where 

^ = z ^ n 2 , E /-w 4Lx(kr, {X}) r, $ . (D.3) 

4V3 A (Kr + A 2 ) d=1 

Here A is the effective gluon mass aimed to incorporate confinement, its value we discuss 
later. The amplitude F ( ^ f _ >x (/c T , {X}) determines the unintegrated gluon density as it is 
introduced in (H), 

J d{X}j2\F^ x (kT,m)\ 2 =4vr^(4,x) . (D.4) 



d=l 



Here 



, = M2 <y- G > ; (D.5) 

MTccG) = !!!i±M + !!i±M + ^±M, (D . 6) 

Oil a 2 O3 

where k\, k 2 , k% and cti, a 2 , 03 are the transverse momenta and fractions of the initial light- 
cone momentum of the projectile gluon carried by the produced c, c and G (see Fig. |T2"| ), 
respectively, and 

k T = ki + k 2 + £3 • (D.7) 

The c-quark spinors £ are defined in @; {X} is the set of variables describing the state X; 
the 10 vertex functions Ti read, 

fi = tf 1 {k x ,a 1 ) Vi{k 2 ^a 2) a z ) ; 

r 2 = -Vi(A;i3, 0:1,0:3) U\{k 2 ,a 2 ) ; 

f 3 = ai Ui_(ki,ai) Vi(k 2 3 - 03^,02,03) ; 

f 4 = -a 2 V 2 (k 13 - a 3 k T ,a 1 ,a 3 )U 2 (k 2 ,a 2 ) ; 

r 5 = 0203^1(^1-^,01)^1(^23,02,03); 

f 6 = -0103^2(^13,01,03) U 2 (k 2 - k T ,a 2 ) ; 
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r 7 = -2a 1 U 1 (k 1 ,a 1 )V 1 {k 23 + a 2 k T ,a 2 ,a 3 ) ; 

f 8 = 2a 2 V 2 (k 13 + axk T , ax, a 3 ) U 2 (k 2 , a 2 ) ; 

f 9 = -2a 2 a 3 U 1 (k 1 - a 1 k T ,a 1 )V 1 (k 23 ,a 2 ,a 3 ) ] 

fio = 2a 1 a 3 V 2 (k 13 ,a 1 ,a 3 )lf 2 (k 2 - a 2 k T ,a 2 ) . 



(D.8) 



Here 



k 13 = a 3 h - aik 3 ; 
^23 = a 3 k 2 - a 2 k 3 ; 
«i + a 2 + a 3 = 1 . 

The matrixes C/j 2 and Vx }2 are defined as follows, 

Ui(ki,ai) = m c a ■ e + (1 — 2ai) e ■ ki + i (e x n) ■ ki ; 

U 2 (k 2 ,a 2 ) = m c (r ■ e + (1 - 2a 2 ) e ■ k 2 + i (e x ft) ■ k 2 ; 
Za6eZ<i4a 

V"i(fc 2 3, "2, "3) = {2a 2 + a 3 ) k 23 ■ e + im c al(n x e) ■ a - ia 3 (k 23 x e) ■ a ■, 

V 2 (k\ 3 ,ai,a 3 ) = (2ai + a 3 ) ki 3 ■ e + i m c a\ {n x e) ■ a — i a 3 (k\ 3 x e) ■ a . 



(D.9) 



(D.10) 
(D.ll) 



(D.12) 



The functions Di in the denominator of ( P-3|) read, 

D 1 = A (A;i, ai) A 2 (k 23 , a 2 , a 3 ) ; 

L> 2 = Ao(k 2 , a 2 ) A 2 (k 13 , a 3 ) ; 

D 3 = A (A;i,Q!i) Ai(fci, k 23 - a 3 k T ,ai,a2,a 3 ) ; 

-D 4 = A (A: 2 , a 2 ) A^, k 13 - a 3 k T , a 2 , a 1: a 3 ) ; 

D 5 = A^kx- k T ,k 23 ,ax,a 2 ,a 3 ) A 2 (k 23 ,a 2 ,a 3 ) ; 

D 6 = Ai(k 2 - k T , k 13 , a 2 , a 1} a 3 ) A 2 (k 13 ,ax,a 3 ) ; 

D 7 = A (kx,ax) Ax(kx,k 23 + a 2 k T ,ax,a 2 ,a 3 ) ; 

D 8 = A (k 2 ,a 2 ) Ax(k 2 ,k 13 + axk T ,a 2 ,ax,a 3 ) ; 

D 9 = A 2 (k 23 ,a 2 , a 3 ) Ax(kx - a x k T ,k 23 ,ax,a 2 , a 3 ) ; 

.D10 = A 2 (fci 3 ,ai,a 3 ) Ax(k 2 - a 2 k T , k 13 ,a 2 ,ax,a 3 ) , 



(D.13) 



where 



A (fc 1 ,ai) = m 2 + k\ — a.\{\ — ax) A 2 ; 
A 2 (fc 13 ,ai,a 3 ) = alm 2 c + a 1 (a 1 + a 3 ) \ 2 + kl 3 ; 
Ai(kx,k 23 ,ai,a 2 ,a 3 ) = (ax — a 3 )(a 3 m 2 + aia 2 A 2 ) + aifc 23 + a 2 a 2 k\ . 
Functions A , A x and A 2 are not independent, but satisfy the relation 



(D.14) 



A (fci,ai) A 2 (£; 23 ,a 2 ,a 3 ) 



ai 



A (kx, ax) Ax(kx, k 23 , ax, a 2 , a 3 ) 



+ a 2 a 3 



A 2 (k 23 , a 2 , a 3 ) Ax(kx, k 23 , a x , a 2 , a 3 ) 



(D.15) 
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The sum of ten amplitudes Eq. ( p.3|) is convenient to split up to four terms, 

10 4 -n 

1=1 i=i v i 



which can be represented using relation (p.!5|) in the following form, 

Mi = ai Ui(ki,oci) fli(h, k 23 ,a 1} a 2 ,a 3 ) + p,x(fa, k 23 - a> 3 k T} a>x,a> 2 ,a 3 ) 
- 2 fix(kx,k 23 + a 2 k T , ai,Qt 2 ,Qt 3 ) 



Mo 



M 3 



-a 2 



foih, k 13 , a 2 , «i, a 3 ) + p, 2 {k 2 , k 13 - a 3 k T , a x , a 2 , a 3 ) 



- 2fi 2 (k 2 ,k 13 + aifcr, ai, a 2 , at 3 ) v 2 {k 2 ,a 2 ) \ 



a 2 a 3 



Ai(/ci, fc 23) ai, a 2 , a 3 ) + Ai(A;i - fcr, &23, "i, "2, 03) 
- 2Ai(fei - aifcr, £23, ai, a 2 , a 3 ) Pi(fc 23 , a 2 , a 3 ) ; 
M 4 = -aia 2 p 2 (A;i 3 ,Q!i,Q; 3 ) A 2 (A; 2 , fc 13 , ax, a 2 , a 3 ) 

+ A 2 (£; 2 - fcr, fci 3 , ai, a 2 , a 3 ) - 2A 2 (£; 2 - ai 2 A; r , k 13 , a x , a 2 , a 3 ) 
The following notations are used here, 



h(k 2 ,a 2 ) 

fai{ki,k 2 3,a 1 ,a 2 ,a 3 ) 
faifa, fa3, ai2, a 3 ) 

Xi(ki, fc 23) ai,a! 2 ,a! 3 ) 
A 2 (A; 2 , ki3, ax, a 2 , a 3 ) 

— * 

pi(k 23 ,a 2 ,a 3 ) 
P2{ki 3 ,a 1 ,a 3 ) 



Ux(kx,ax) 
A (kx,ax) 
U 2 {k 2 ,a 2 ) 
A (k 2 ,a 2 ) 



Vx(k 2 3 ) a 2 ,a3) _ 
Ai(/ci,A; 2 3,ai,a 2 ,a3) 

Ax(k 2 ,kx3,ax, a 2 ,a 3 ) 



^l(fcl,Q!l) 

Ai(/ci,fc 23 ,ai,a 2 ,a 3 ) 

?7 2 (fc 2 ,Qi 2 ) 

Ai(A; 2 ,A;i3, ai,a 2 ,a 3 ) 



Vl(fc23,Q2,Q3) 

A 2 (/c 23 ,a 2 ,a 3 ) 

V2(fcl3,Q!i,Q3) 
A 2 (fcl3, OCX, Oi 3 ) 
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Apparently, the amplitude operators in Eqs. ( D.17 )- (|D.20| ) vanish, — > when %t — > 
what guarantees infra-red stability of the cross section of gluon radiation. 

Since the cross section steeply falls down at large k\ and Mi{k T = 0) = we can expand 
the amplitudes at small kx, 



M t {k T ) 

Using Eqs. (|DT^ )- (|LT20l ) we find, 



_d_ 

dki 



Mi{kj 



(D.25) 



J k T =0 



8M l 
dk T 


k T =0 


~ — «i V\ 


dM 2 
dkx 


k T =0 


dpL 2 

7H at -=r~ 

dk 13 


dM 3 
dkx 


k T =0 


« -a 2 Q!3 


<9M 4 
dkx 


k T =0 





(a 3 + 2a 2 ) 



dk 23 

■ z> 2 (a 3 + 2 a x ] 



gAi 

<%2 



Pi(l -2ai; 
l-2a 2 ) . 



(D.26) 
(D.27) 
(D.28) 
(D.29) 



Using the obvious relations 



1 - 2a t 
l-2a 2 



a 2 — a i + a 3 
a 1 — a 2 + a 3 



(D.30) 



we conclude that the terms M 3 and M4 in the amplitude of cc pair production are negligible 
compared to M\ and M 2 - Indeed, since the cc must be projected on the wave function of 
the heavy charmonium, one can make use of the fact that it is a nonrelativistic system, and 
that the radiated gluon is predominantly soft, i.e. 



\k\ - k 2 \ < m c , 
|«i — a 2 \ <C 1 , 
a 3 « 1 ■ 



(D.31) 



These conditions also help to simplify essentially the Eqs. ( p.26|) and ( p.29|) for M\ and 
M 2 . First of all, for vanishing a 3 one can neglect the terms containing the Pauli matrixes 
in the definition (|D.12| ) of V\ and V 2 which now commutate with U\ and U 2 . 

Further, for small a 3 Eq. (|D.9| ) leads to k\ 3 ~ —a\ k 3 , k 23 « —0:2 fc 3 . Since from ( |D.31| ) it 
follows that «i ~ a 2 ~ 1/2 eventually we arrive at &i 3 ~ ^23 ~ ~k 3 /2. Then the propagator 
A x can be represented as 



Ax 

T 2 

M 



(kl+T^+OMh-k,)] 



X 2 + a 3 M 2 



2m c 



(D.32) 
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The approximations done we arrive to a simplified form of the amplitude, 



M 



m c a ■ e + i (n x k\) ■ e m c a ■ ei {n x k 2 ) ■ e 



x 



m 2 + k\ 
2(k 3 -k T )-e 2k 3 -e 



ml + k 2 



(D.33) 



Here we suppressed in the propagator the term A 2 /4 ~ Aq CD /4 which is tiny compared to 
Using Eq. ( |D.33| ) the amplitude ( |L).3| ) takes the form 

A^(ki2, k 3 , k T , a 3 ) = d 2 bd 2 sd 2 r A^ (b, s, f, a 3 ) exp i(b ■ k T + s- k 3 + r- k 12 , (D.34) 



m 



where 



A%(b,s,r,a 3 ) 



V3 



J2f M ^(r,a = l/2) 



<i=i 



7 



(b + s,{X})-^(b- r -,{X} 



- $ 



qG 



7 W (& + S -,{X})- 7 ( d )(6+^{X}) 



(D.35) 



Here (r*, a) is the distribution amplitude for a cc fluctuation of a gluon given by Eq. (|7|), 
and Qgcis) is the distribution amplitude for a quark-gluon fluctuation of a quark which has 
the form, 

^c(r) = - x Ue-V r K (rr) . (D.36) 

7T V O 

The profile function 7^(6, {X}) in Eq. (|D]35|) is related to the amplitude F^ d \k T , {X}) 
by Fourier transform, 

d 2 k T 



7 W (6, {X}) 



GN^X 



(k T , {X}) . 



(D.37) 



2tiV6-' /# + A 2 

This profile function is related to the unintegrated gluon density T{k,T-, x) and the dipole 
cross section <Tq q (r, x) by the relation, 



J d 2 bdv x j2[i {d \b + f,{x})- 1 {d \b,{x}) 



Air 



a. 



d 2 kn 



kl + X 2 



1-e 



T{k T ,x) = cr qq (r,x 



(D.38) 



Now we can calculate the amplitude of x production accompanied by gluon radiation 
which corresponds to the graphs depicted in Fig. [12] with the cc pair projected to the wave 
function of x, 



ii/j 



A* b (k T ,k 3 ) =J2 I dad 2 bd 2 sd 2 r$^(r,a) A^(b,s,r) exp \i{b ■ k T + s ■ k 3 ) 
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(D.39) 



Note that the typical value of r ~ l/m c is much smaller than the mean separation 
b ~ 1/Aqcd- Therefore, we can use the approximation 7^(6 ± r, {X}) ft 7^(6, {X}) in 
Eq. (|D.35| ), then the amplitude takes the form which is known for the radiation of a 



heavy photon by a quark, i.e. the amplitude of Drell-Yan process [[26 



A* b (k T ,k 3 ,a 3 ) ft ^ f d 2 bd 2 s^(s,a 3 ) exp i(b ■ k T + s- k 3 ) 

1 mm 

7 (d) (6 + s, {X}) - 7 (d) (&, {X})] , (D.40) 



where \I/(s, 03) defined in (|50l) is the effective distribution amplitude for s x ~ G fluctuation 
of a gluon, which is the analog for the 7* q fluctuation of a quark. 

The two terms in square brackets in (|D.40| ) correspond to the two diagrams in Fig. |6] 



which are the same as for Drell-Yan reaction if to replace the virtual photon by x an d the 
incoming and recoil quarks by gluons. Therefore, it is natural to expect the corresponding 
cross section to have the same factorized form as for Drell-Yan reaction. Indeed, squaring 
expression Eq. ( p.40| ) and using Eq. ( |D.38| ) we arrive at the expression Eq. ([49]) for the cross 
section. 
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